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Optimal pair density functional for description of nuclei with large neutron excess 
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Toward a universal description of pairing properties in nuclei far from stability, we extend the 
energy density functional by enriching the isovector density dependence in the particle-particle chan- 
nel (pair density functional, pair-DF). We emphasize the necessity of both the linear and quadratic 
isovector density terms. The parameters are optimized by the Hartree-Fock-Bogoliubov calculation 
for 156 nuclei of the mass number A — 118 — 196 and the asymmetry parameter (A'' — Z)/A < 0.25. 
We clarify that the pair-DF should include the isovector density dependence in order to take into 
account the effect of the isoscalar and isovector effective masses in the particle-hole channel con- 
sistently. The different Skyrme forces can give the small difference in the pairing gaps toward the 
neutron drip line, if the optimal pair-DF consistent with the particle-hole channel is employed. 

PACS numbers: 21.10.Re, 21.60.Ev, 21.60.Jz 



I. INTRODUCTION 

The energy density functional (EDF) theory provides 
a comprehensive microscopic framework for description 
of bulk nuclear properties, low-lying excitations, giant 
vibrations, and rotational excitations [l|. From the pi- 
oneering work by Vautherin and Brink diverse en- 
deavors have been made for finding the best EDF aiming 
at the description of the nuclear properties across the 
mass table. For example, the Skyrme functional for the 
particle-hole (p-h) channel has been improved by tak- 
ing into account the incompressibility modulus of nuclear 
matter 3], the spin- and spin-isospin channels |4], the de- 
formation properties Q , the spin-orbit terms 6] , and the 
isospin properties [t^. Efforts to include the new terms 
such as the tensor terms are also being made (for the 
recent situation, see Ref. Q). 

The particle-particle (p-p) channel of the EDF (pair 
density functional, pair-DF) is also an indispensable cl- 
ement for description of nuclear systems Q. The study 
of the nuclear matter predicts a very weak ^Sq pairing 
at the normal density, and the pairing correlation in fi- 
nite nuclei is considered to be nuclear surface effects [l^] ■ 
The induced pairing interaction due to phonon exchange 
also enhances the surface effect Ell, [13 ■ These facts 
suggest the density dependence of the effective pairing 
force. 

The standard parametrization of the pair-DF has the 
isoscalar density {p ~ p„ + pp) dependence only [3, [3, 
[Ts. . .16. , ,17] . The coupling constant should be constrained 
by the requirement to reproduce the experimental data 
such as masses, low-lying excited states, and rotational 
properties. However, the functional form of the density 
dependence is still under discussion [l^ . 

In nuclei near the /3-stability line, the effect of the p-h 
field characterized by the Fermi energy is much stronger 
than the p-p field. Therefore the pairing correlations 
can be treated within the BCS approximation 0, [l9| . 



On the other hand, the strengths for the p-p and p-h 
channels become cornparable in magnitude for weakly- 
bound nuclei [1, [H, Hfl, HH [H, [H, [2i| ■ Therefore it is 
desirable to constrain the functional form of the pair-DF 
by using the experimental data of unstable nuclei. 

The isovector density (pi — pn — Pp) dependence can 
have sizable effects in nuclei apart from the /3-stability 
line. In Ref. psj . the linear pi terms were introduced so 
as to simulate the neutron pairing gaps in symmetric and 
neutron matters obtained with either the bare interaction 
or the interaction screened by the medium polarization 
effects. It was pointed out that the pairing properties 
in semi-magic nuclei can be better described by the pi- 
dependent pair-DF than that without pi terms 26|. 

We also recognized the importance of the linear pi term 
in the pair-DF [24i|. By performing the Hartree-Fock- 
Bogoliubov (HFB) calculation with various coupling con- 
stants of the p and pi terms, we emphasized the strong 
sensitivity to the pairing properties and the influence on 
rotational excitations in deformed nuclei near the neu- 
tron drip line. 

In principle, it is desirable to derive the pair-DF from 
the bare interaction based on the microscopic pairing the- 
ory including both the medium polarization effect and 
the surface phonon coupling effect in finite nuclei. How- 
ever, it seems very difficult at present in spite of recent 
progress toward this direction [H, [13, [H, [21, [H [Ml, [H, 
|3l,[3i,[3i]. 

In this paper, we extend the pair-DF by including 
the linear and quadratic pi terms based on the phe- 
nomenological considerations. The pair-DF is designed 
so as to reproduce the dependence of pairing gaps on 
both the mass number A and the asymmetry parameter 
a = (N — Z)/A. The parameters in the pair-DF are op- 
timized so as to minimize the root-mean-square (r.m.s.) 
deviation between the experimental and calculated pair- 
ing gaps. The necessity of the pi dependence in pair-DF 
is emphasized in connection with the the effective mass 
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FIG. 1: (Color online) Experimental neutron pairing gaps 
(upper panel) and those of proton (lower panel) in the region 
of N, Z > 40 (except for nuclei with either Z — 50, 82, or 
A*' — 50, 82, 126) are shown as a function of a. The pairing 
gaps are divided by A^*'. The error bar represents the r.m.s. 
deviation from the average trend A^^^\a) for each a with 
Aa = 0.02 interval. 
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We extend the analysis with up-to-date measured 
masses in the wider mass region of TV, Z > 40 (except 
for nuclei with either Z = 50,82 or TV = 50,82,126) 
[sj • The result is shown in Fig. [T] The average A- and 
a-dependence is determined for the neutron and proton 
pairing gaps separately by x^-fitting; 

^(oxp) 



A (cxp) 
P 



(a)/Af) 



^(0) , 



1 - 7.74^2 



p(O) 



= 1 - 8.25 



(1) 



(2) 



with a1^^ = 6.75/A1/3 MeV and A^"^' = 6.36/^1/3 MeV. 
Here the experimental pairing gaps are extracted by the 
odd-even mass difference with the three-point staggering 
parameters 18]. 

The Coulomb force is an important building block of 
nuclear systems. The 20 - 30 % reduction of Ap by the 
self-consistent treatment of the Coulomb force was re- 
ported in Ref. jsl] . The authors of Ref. [s^ also arrived 
at the same conclusion by performing the HFB calcula- 
tion with the non-empirical pair-DF. On the other hand, 
the experimental evidence is unclear in our analysis. The 
ratio is A^""P^(a)/Ai""P)(a) « 0.94(1 -0.51a2) > 0.91 for 
< a < 0.25. This is smaller than the uncertainty of our 
analysis about 10 % shown by error bars in Fig. [TJ The 
elaborate investigation is required to clarify the Coulomb 
effect. Therefore we neglect this effect in our analysis and 
leave it as an open problem in the future study. 



parameters. 

This paper is organized as follows. In Sec|lTl we briefly 
review the global properties of pairing gaps. In Sec lllll 
our pair-DF including the linear and quadratic pi terms 
is introduced. In Sec. IIVI we point out the drawback of 
the pair-DF without the pi terms. In Sec. IVland lVIi we 
investigate the role of the pi terms. The parameters in 
the pair-DF are determined by the HFB calculation for 
156 nuclei of A = 118 - 196 and a < 0.25. We clarify the 
close connection between the pair-DF and the effective 
masses by the extensive analysis with 13 Skyrme param- 
eters. The choice of the pairing strength is discussed in 
Sec. I VIII The conclusion is drawn in Sec. IVIIII 



II. GLOBAL TREND OF PAIRING GAPS 

We construct the pair-DF so as to reproduce the A- 
and a-dependence of pairing gaps. In Ref. [11], Vogel et 
al. pointed out that the experimental pairing gaps in the 
region of 50 < Z < 82 and 82 < TV < 126 can be well 
parametrized by A (a) = (l — 6.1 a^) A^^^ Here the A 
dependent part is given by A^^^^ = 7.2/^^/^ MeV. This 
average ^-and a-dependence holds for both neutron and 
proton pairing gaps. 



III. MODEL 

A. Parametrization of pair-DF 

We extend the pair-DF by including the linear and 
quadratic pi terms in the following form; 



pair \ ' } — ^ 



T—n,p 



9r [p,Pi] {Pt {r)Y 



with 

9t [p,Pi] = 1 



p{r) T3pi{r) 
Vo Vi 



Pi jr) 
Po 



(3) 



■(4) 



Po Po 

Here t = n (neutron) or p (proton), and po = 0.16 fm~'^ 
is the saturation density of symmetric nuclear matter. 
The T3 = 1 (n) or —1 (p) in the linear pi term is in- 
troduced so as to preserves the charge symmetry of the 
pair-DF. In nuclei with large a, the pi terms produce 
two effect for pairing correlations. The one is the volume 
effect inside the nucleus, which is relevant to all nuclei. 
The other is the skin effect in nuclei apart from the (3- 
stability line. 

The pair-DF with 770 — 0.5 and 771 = 772 = is one of 
the current standard parameterizations called the mixed- 
type pairing force. This pairing force reproduces the av- 
erage A dependence of pairing gaps [l^. We also justify 
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this choice in Sec. |Vl Therefore, we fix rjo = 0.5 unless 
otherwise noted. 



B. Setup 

We use the standard Skyrme interaction for the p-h 
channel in the HFB calculation. The Skyrme SLy4 [7:] 
parametrization is mainly used. In Sec. IVIl we will ex- 
tend our analysis with 13 Skyrme parameters. 

For the determination of rji and 772, we perform the 
Skyrme-HFB calculation for 156 ground states of even- 
even, open-shell nuclei in the region oi Z = 56 — 76, and 
either iV = 56 — 76 or 88 — 120, which covers the range 
of < a < 0.25. We utilize the computer code of the 
Skyrme-HFB calculation developed by M. Stoitsov et al. 
[4q |. Starting from the spherical, prolate and oblate ini- 
tial conditions, the lowest energy solution is searched in 
the space of axially symmetric quadrupole deformation. 

We estimate the r.m.s. deviations between the experi- 
mental and calculated pairing gaps in order to optimize 
the 771 and 772. The neutron and proton r.m.s. deviations 
are defined by 
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The total r.m.s. deviation is also given by 

1/2 
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Here Nr'^^^'' is the number of existing data of A^"'^'''' 



in the region of the present investigation; N, 
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(cxp) 



and Nj 
by |4il,|4i,|43 
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84. The theoretical pairing gap is defined 



drprir)hrir)/ / drpr(r). 



A, = - 

when the local pairing potential is given by 
d 



hr{r) 



dpr{r) 



We extract the coefficients Cr which represent the pairing gaps are A- 
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Vo[A„(^5«Dy)] 


""tot 0"p 


r'(i) r-^") r*-^^ 


50 
75 


-346.5 
-320.0 


0.17 0.16 0.17 
0.17 0.16 0.18 


1.08 9.42 1.00 8.13 
1.07 9.26 1.01 8.44 



is constrained by the aL""^' of 



TABLE I: The cutoff quasiparticle energy Ecut dependence 
of the r.m.s. deviations [MeV] and the coefficients C|*' are 
listed. The parameters {rjo , rji , 1^2) = (0.5,0.2,2.5) are fixed 
The strength Vo [MeV fm 
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FIG. 2: (Color online) Neutron and proton pairing gaps ob- 
tained with 770 — 0.5 and ?7i = 772 = are plotted as a function 
of a. The pairing gaps are divided by Ai^\ The eoNd and 
70 Yb isotopes possessing the large proton shell gaps are indi- 
cated by circles with horizontal and vertical bars respectively 
in the bottom panel. 



(exp) 



average a-dependence of pairing gaps, 

A.(a) = (c(°)-C«a2)A(^), (9) 



1.17 MeV and A 



(cxp) 



= 0.98 



MeV, which are close to Ai°"P^(a) 



by x^-fitting analysis for Ar of the 156 nuclei. Here Ar^-* 
is the same quantity determined for Eqs.iJT]) and ((2]). 
For each set of (770, ??2), the strength Vq is fixed so 

as to reproduce the A^°'^^'' of ^^^Dy. We use the abbre- 
viation Vo[A„(^^^Dy)] for this choice. This nucleus has 
quadrupole deformation /3 ~ 0.28 The experimental 



1.04 MeV and 

A^^'^P^ (a) = 0.96 MeV estimated by Eqs.HI) and ©. The 
justification of Vq will be discussed in Sec. IVIII 

The cutoff quasiparticle energy Ecut = 50 MeV is fixed 
in this paper. We checked the dependence of ctt and Ct 
on the i^cut in Table |T| The results with Ecut — 50 and 
75 MeV agree within a few percent accuracy. Here the 
parameters of the pair DF are fixed to be (770,771,772) = 
(0.5,0.2,2.5), which are the optimal choice (See Sec. |V| . 
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FIG. 3: (Color online) The coefficient cl"' obtained with r\\ = 
r;2 = is shown as a function of r;o. The result with 772 = 2.5 
and ?7i(r;o,7?2) is compared. Here 171(770, 772) is the value of r\\ 
minimizing (Jtot for each (770,772) with Vo[A„("'^^'^Dy)]. 
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FIG. 4: (Color online) The same with Fig. [3] but for the co- 
efficient C^^'. 



IV. PAIR-DF WITHOUT pi TERM 

We show the drawback of the pair-DF without the 
p\ terms. The pairing gaps obtained with r/o = 0.5, 
?7i = ?72 = 0, and the strength Vb[A„(i56Dy)] = -324.0 
MeV fm-3 g^j.g plotted in Fig. H The extracted cf^ 
and ci^'' are shown by the dashed lines in Figs. [3] and U) 
The A„ and Ap are almost a-independent. The coeffi- 
cient ci^"* : 

value Ci^oxp 
^(1) 



1.11 is much smaller than the experimental 
7.74. Although the cf' = 3.74 is larger 

than Cn' = 1.11, this is due to the collapse of the pair- 
ing gap in weak pairing region. Actually, it would be 



(1) 



1.38. 



if we restrict the data to Ap > 0.25 MeV. 



Here the eoNd and 70 Yb isotopes possessing the large 
proton shell gaps are indicated by circles with horizon- 
tal and vertical bars respectively in the bottom panel of 
Fig. [21 The coUagse of Ap is a drawback of the mean-field 
approximation [45|, and can be overcome by performing 
the particle number projection (PNP). The improvement 




FIG. 5: (Color online) The same with Fig. [3]but for a^. 
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FIG. 6: (Color online) The r.m.s. deviations are shown as 
a function of 771. The results with either 772 = or 2.5 are 
compared. 

for the Cr^-* values, however, can not be expected by the 
PNP procedure [46], because the PNP method does not 
have any specific isovector effect. Therefore, we neglect 
the effect of PNP in this study. 
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Cp'' = 0.67 is smaller than Cn' due to the quenching of 
Ap attributed to the neutron skin effect [24 1: The neu- 
tron skin reduces the overlap between the form factor 
[1 - %p('')/Po] and /5p(r) in Eq.®. 

The quenching of Ap due to the neutron skin effect 
becomes stronger with larger 770 [^Jj. The a-r is shown 
as a function of 770 by the dashed line in Fig. [5l Because 
the cTp rapidly increases with 770, the minimum of dtot is 
absent. In addition, the cf^ and Cr^"* remain small if 
restricted to 771 = ?/2 = (See Figs. [3] and H]). 



V. ROLE OF pi DEPENDENCE 

It is possible to compensate the quenching of Ap by us- 
ing a str onger p airing strength for proton (for example, 
Ref. fisl-ISSllSq). However it violates the charge symme- 
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FIG. 7: (Color online) The r.m.s deviations at rii = 0.2 are 
shown as a function of 772- 




FIG. 8: (Color online) The coefficients at 7?i = 0.2 are 
shown as a function of Tj2. 
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FIG. 9: (Color online) The same as Fig. [2] but with 
(77o,?7i,^2) = (0.5,0.2,2.5). 
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-346.5 
-322.0 


0.17 0.16 0.17 
0.34 0.30 0.37 


1.08 9.42 1.00 8.13 
0.76 7.87 0.69 6.11 


Exp 






1.00 7.74 1.00 8.25 



TABLE II: The r.m.s. deviations and the coefficients Ct ^ 
obtained with the SLy4 force and the optimal parameters 
(770 , ?7i , ) ~ (0.5,0.2,2.5) are shown. The results with 
yopt(def), l/opt(sph), yo[A„(i56Dy)], and V^o[A„(i20Sn)] are 



compared. The experimental values of C|*' are also listed 



try of the pair-DF. This is the important symmetry in 
the theoretical framework, and indispensable for global 
description of pairing properties from neutron to proton 
drip line. 

This consideration leads to introduction of the linear 
Pi term in Eq. This pair-DF preserves the charge 
symmetry. The pi term induces the difference of the 
neutron and proton pairing strengths automatically [2^ . 
The Up has the minimum value at rji = 0.15, while the (T„ 
is almost constant as a function of 771. This is shown by 
the dashed hues in Fig.[6l We obtain the Cn^ = 0.93 and 
Cp"-* — 0.83 with r/i = 0.15 and 772 = 0, which are better 
than those with ?7i = 772 = 0. However, the Cn ^ = 3.75 
and Cp^^ = 1.89 at r/i — 0.15 do not improve as a function 



of 771. 

The quadratic pi term in the pair-DF improves the 
r.m.s. deviations. To see this, the r.m.s. deviations are 
plotted as a function of 7/1 while keeping 772 = 2.5 in 
Fig. [HI Those with rji = 0.2 are plotted as a function of 
772 in Fig. [T] The parameter set of (771,772) = (0.2,2.5) 
simultaneously gives the minimum values of crtot , and 
CTp. They are smaller than those at 772 = 

The quadratic pi term also improves the a- 
dependence. The Cr"'' and cj^"* at 771 = 0.2 are plot- 
ted as a function of 772 in Fig. [8] The ci^^ stays 
around 1.0, while the Ct^^ increases linearly and reaches 

Ct ' ~ Cr,oxp ~ 8 at 772 = 2.5. 

The pairing gaps obtained with (771,772) = (0.2,2.5) 
are shown in Fig. [5] The r.m.s. deviations and the coef- 
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FIG. 10: (Color online) The value of ?7i minimizing the a^ot 
for each (»7o,??2) with Vb[A„(^^®Dy)] is plotted as a function 
of Tjo- The 7)2 = 2.5 is fixed. 



ficients C'^^ are listed in Table HTl We see the significant 
improvement compared to those with rji = r]2 ~ 0. 

The optimized set of (771,772) gives the justification of 
the mixed type pairing force (770 = 0.5). The ar with 
772 = 2.5 and rii{riQ,ri2) is shown as a function of 770 by 
the solid line in Fig. [5] The improvement over the choice 
Vi = V2 = Oj especially the large reduction of dp, is 
obvious. Therefore, the minimum of crtot can appear at 
7/0 ~ 0.5. Here 771 (7/0,772) is the value of 771 minimizing 
cTtot for each (770,772) with Vb[A„(i56Dy)]. The 7/1(770,772) 
at 772 = 2.5 is shown as a function of 770 in Fig. 1101 



The coefficients C)- 
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and Ci^''> 



with 772 = 2.5 and 



771 (770 , 772 ) are shown as a function of 770 by the solid lines 
in Fig. [3] and [H The Ct"'' is insensitive to 770 , while the 
Ct^^ becomes close to the experimental value at 770 0.5. 
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FIG. 11: (Color online) The strength Vb reproducing the 
A^'^'^''^ of ^''''Dy for each Skyrme force is plotted in relation to 
m*/m. The results with 772 = 2.5 and the optimal r/i in Table 
mil are compared to those with r/i = 772 = 0. The 770 = 0.5 is 
fixed. 
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VI. EFFECTIVE MASS AND pi-DEPENDENCE 
OF PAIR-DF 



FIG. 12: (Color online) The optimal value of 771 minimizing 
(Jtot with (770,772) = (0.5,2.5) and Vo[A„(^^^Dy)] is shown in 
relation to m/m*. See text for details. 



A. Isoscalar and isovector effective masses 

Pairing correlations are sensitive to the single-particle 
structure around the Fermi level. For a suggestive ex- 
ample, the pairing gap is a function of gG and given by 
A cx e~^/^*^ for gG ^ 1 in the schematic model of the 
seniority pairing force with the strength G and the uni- 
form single-particle level density g ^M,- On the other 
hand, the effective mass has a strong infiuence on the 
single-particle energies. The average level density is pro- 
portional to the effective mass [l3|- Therefore, we expect 
the close connection between the effective mass and the 
pair-DF in order to reproduce the global trend of the 
experimental pairing gaps. 

For the investigation, we extend our analysis with 13 
Skyrme parameterizations; SkM* d], SGII [|, LNS [13, 
SkP BSkl7 0, SkT6 [ipL SLy4, SLy5 [7], Skll, 
Skl3, SkM [el], SkO, and SkoTHg. 



The effective mass of the Skyrme force is given by 



bip - h\pr 



2771* (r) 2777 

Zm 777* 



(10) 
(11) 



with /(r) = Pi/p and A777i(r) = ra/m* — m/m* pi [s^. 
The isoscalar and isovector effective masses are defined 

by 



m 



mj(r) 

m 
mlir) 



= 1 



2777 / b[ 
27^7 



hip — 1 ~\- K. 



(12) 
(13) 



The 777* is directly connected to the enhancement factor 
K of the Thomas- Reiche-Kuhn sum rule ^jj. We esti- 
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SGII 


0.670 0.786 -0.219 
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0.325 


-321.3 


0.15 0.15 0.16 


1.09 9.38 


0.95 8.10 


LNS 


0.727 0.825 -0.164 


1 


0.325 


-322.2 


0.17 0.13 0.21 


1.09 11.58 


0.98 8.95 


SkP 


0.732 1.000 -0.366 


1 1 


0.300 


-268.0 


0.16 0.19 0.12 


1.07 9.05 


0.92 8.24 


BSkl7 


0.780 0.798 -0.028 


1 1 


0.200 


-313.5 
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1.08 8.68 
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SkI4 


0.800 0.649 0.290 


-0.985 


0.275 


-364.5 


0.18 0.16 0.20 


1.02 8.25 


0.97 7.57 


SkO 


0.851 0.896 -0.058 


-1.125 


0.250 


-290.5 


0.17 0.18 0.17 


1.05 6.91 


0.91 6.24 


SkO' 


0.871 0.896 -0.032 


-0.576 1 


0.200 


-289.0 


0.16 0.16 0.15 


1.03 6.62 


0.89 5.64 


SkI3 


0.800 0.574 0.493 





0.125 


-397.4 


0.20 0.20 0.20 


1.04 7.26 


0.98 8.63 


SkT6 


1.000 1.000 0.000 


1 1 


0.075 


-268.9 


0.17 0.21 0.11 


1.11 7.93 


0.92 6.97 



TABLE III: The parameter set of the optimal pair-DF for each Skyrme parameterization is listed. The optimal value of 771 
minimizing atot with the strength Vb[An(^^^Dy)] [MeV fm~^], the r.m.s. deviations [MeV], and the coefficients Ct'' are shown. 
The parameters (rjo, rj2) = (0.5, 2.5) are fixed for them. The effective masses m* and m* at the saturation density of symmetric 
nuclear matter, the difference Ami ~ m/m* — m/m^, and the Wq/Wq and rjj of the spin-orbit potential are also listed. 
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FIG. 13; (Color online) The same with Fig. 1 121 but in relation 
to Ami. 



mate the m* and m* at the saturation density of the 
symmetric nuclear matter. The m* is a key parameter 
which determines the splitting of the neutron and proton 
effective masses as a function of /. The parameters bi 
and b[ are given by 

bi = ^[ti{2 + xi)+t2{2 + X2)] (14) 



- [ii (1 + 2x1) - <2 (1 + 2x2)] . (15) 



The spin-orbit potential W,- (r) also has the pi depen- 
dence. The WT-(r) of Skyrme DF is defined by 



W,(r) = -(WoVp + W-^VpO + ^jW;? 

= hr + -rW/^ + ^3-rVpi + 77jW;^ 



where W;?(r) = CqJ + rsC/Ji with the parameter rjj 
of either or 1. The J (Ji) is the isoscalar (isovector) 
spin-current density. The Cq and C/ are given by ii, t2, 
xi and X2 of the Skyrme parameter [HI]. Most Skyrme 
functional have the spin-orbit terms with Wq = Wq. 
However, the SkI4, SkO, and SkO' have the generaUzed 
pi dependence by introducing the different strengths Wq 
and W^. The SkI3 has = 0. 

We search the optimal value of rji which minimizes dtot 
under the conditions; 1) the fixed (r/o,?72) = (0.5,2.5), 
and 2) the strength Vq reproducing the An^'^^ of -'^^^Dy. 
The numerical uncertainty is drji = 0.025. The results 
are summarized in Table [Till 

The strengths Vb reproducing the An^^^ of ^^^Dy are 
plotted in relation to m*lm in Fig. 1111 For 771 = 772 = 0, 
the Vq increases linearly. This trend is in agree with the 
general consideration that the pairing strength should 
be increased if the level density is low. The relation is 
given by Vb = -505.05 + 2MAlm*Jm MeV fm'^ with 
the correlation coefficient r — 0.99, except for SkT6 and 
SkI3 (See AppendixIXlfor the procedure of the correlation 
analysis). In general, the linear correlation disappears 
with the pi terms due to the 771 and 772 dependence of 
Vq. However, it is interesting to mention that the linear 
correlation is recovered with 772 = 2.5 and the optimal 771. 
The extracted correlation is Vb = —531.45 + 266. 4677i*/?7i 
MeV fm-3 with r = 0.99. 

The optimal values of 771 are shown in relation to m/rrv^ 
in Fig. 1121 The linear correlation between 771 and m/m'^ 
is obvious, irrespective of the choice of Wq and rjj. The 
extracted relation is given by 



7^1 = -0.340 + 0.464 



m 



or, in terms of the enhancement factor 
7^1 = 0.124 + 0.464 k, 



(16) 



(17) 
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except for SkI3 and SkT6. The correlation coefficient 
between 771 and m/m* is r = 0.85. This indicates that 
these parameters is almost linearly dependent. The pos- 
sible reason for the deviation of SkI3 and SkT6 is the 
special assumption on the Skyrmc DF. The SkT6 sets 
TO* — m* = TO by definition. The SkI3 neglects the pi 
term in the spin-orbit potential by setting Wq — 0. 

During the optimization of the pair-DF, the 772 = 2.5 
for SLy4 is used for other Skyrme parameters to avoid 
the huge computational task. However, the improvement 
obtained by the optimization of the parameter r]2 should 
be small. The effect can be estimated as follows: We 
define the r.m.s. deviation of C'^' by 



1.0 



{a 



(1) 



(18) 



Here ci]c^p is the experimental value of Eqs. ^ and 
(HI). The AC(i) = 1.5 is obtained by taking the average 
<> over T = n,p and the 13 Skyrme parameters. If the 
linearity A772 « AC^^y2.3 of Fig. [8] and the parabolic 
approximation for dtot as a function of 772 from Fig. [7] are 
assumed for the other Skyrme parameters, the expected 
improvement for ctot is about 0.002 MeV. In addition, 
we will show that the difference in the pairing gaps for 
different Skyrme forces is small if the pair-DF with 772 = 
2.5 and the optimal rji is used in Sec. IVIII 

From the present analysis, we conclude that the pair- 
DF should include the pi dependence in order to take 
into account the effect of the 771* and to* for the global 
description of pairing correlations. 



B. Attii dependence 

The effective masses 771* and to* have strong correla- 
tion with the Vq and 771 of the pair DF respectively. On 
the other hand, the splitting of neutron and proton ef- 
fective masses directly depends on the local asymmetry 
parameter /(r) through the combination of to* and 771*; 
A7771 = m/m* — m/ml- If A77t.i is negative, the m* 
(to* ) is a decreasing (increasing) function of /. In order 
to compensate the effect, the larger 771 is necessary so 
as reproduce the same magnitude of the A'^^^\a) and 
Ai°''P^(a) (See Eqs. and ©). On the other hand, 
the smaller 771 is required for the positive A7ni. 

This correlation can be seen in Fig. [T31 The extracted 
correlation is 



771 = 0.261 - 0.193A7711 



(19) 



with r = —0.63, except for SkT6 and SkI3. Although this 
correlation is weaker than that between 771 and m/m* 
due to the scattering of the m*/m value, it is meaningful 
to conclude the linear dependence between 771 and A7711. 
The 771 « 0.26 at Atoi k, can be the rough estimation 
of the 771 compensating the artificial Ap suppression due 
to the neutron skin effect. 
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FIG. 14: (Color online) The r.m.s. deviations with SLy4 
force and (770, 7?i, 772) = (0.5,0.2,2.5) are shown as a function 
of Vo. The strengths K>pt(def), K>pt(sph), l/o[A„(i56Dy)], and 
Vb[A„(^^''Sn)] are indicated by arrows. 



K„p,(def) ---■Fo(A„('''Dy)) 

- - - F„p,(sph) ro(A„("°Sn)) 
* exp. 



Z=82 




0.2 0.3 
a={N-Z)lA 
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FIG. 15: (Color online) The neutron pairing gap of Sn and 
Pb isotopes with SLy4 force and (770,771,772) = (0.5,0.2,2.5) 
are shown as a function of a. The results with the strengths; 
Vopt(def), l/opt(sph), yo[A„(i5''Dy)], and Vo[A„(i20Sn)] are 

compared. The experimental trend An^^^ (a) is shown up to 
a < 0.25 together with the experimental data. 



VII. CHOICE OF Vo 

It is desirable to optimize the strength Vq for exper- 
imental data in wide region of nuclear chart. How- 
ever, the procedure demands heavy computational ef- 
forts. Therefore the Vq is usually fixed so as to repro- 
duce a pairing gap of specific nucleus. Several authors 
adopted the aI'^'^p^ of ^^OSn [H [H, The strength 
I/o[A„(i20gj^)] ^ -322.0 MeV fm'^ with SLy4 force and 
(77^,772) = (0.2,2.5) gives CTtot = 0.34 MeV (see TablcHH. 
On the other hand, if the Vq is fixed in deformed re- 
gion, for example, yo[A„(i^^Dy)] = -346.5 MeV fm'^ 
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* exp. 
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FIG. 16: (Color online) The same with Fig. [15] but for the 
proton pairing gaps of A'^ = 50 and 82 isotones. 
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FIG. 17: (Color online) The difference in the pairing gaps for 
the different Skyrme parameters. See text for details. 



the r.m.s. deviation reduces to ctot — 0.17 MeV. The Vb 
dependence of the Ctot is shown in Fig. [TH The optimal 
value yopt(def) = -344.0 MeV fm"^ for the deformed 
nuclei is close to the Vo[A„(i5SDy)]. 

In Figs, [m and [TH the A„ of Sn and Pb isotopes 
50 and 82 isotones obtained with 
,772) = (0.2,2.5) are shown. The re- 



and the Ap of 
SLy4 force and (771 

suits with yopt(def), yo[A„(i56Dy)], yo[A„(i20Sn)] are 
compared. The difference of the results with yopt(def) 
and Vb[A„(^^^Dy)] is negligible along the isotopic and 
isotonic chains. However, the choice of Vn[An("'"^^Dy)] 
overestimates the experimental pairing gaps [37| in these 
spherical nuclei. The strength optimized only for the 
spherical nuclei is yopt(sph) = -308.0 MeV fm"^. This 



is 10.4 % weaker than T4pt(def). 

It is an open problem to construct the pair-DF which 
allows us to describe the pairing properties along the 
chains of semi-magic nuclei at the same quality achieved 
for deformed region [El, [H, [H^ • The authors of Ref. (5^1 
considered that the overestimation in spherical nuclei 
may be partly attributed to the elfect of the particle num- 
ber fluctuation. They showed that the HFB calculation 
with the approximate particle number projection using 
the Lipkin-Nogami method improves the agreement with 
experiment for spherical nuclei. We do not discuss this 
point further in detail, and the choice of Vb[A„(^^^Dy)] 
is employed in this work. 

The strengths Vq for other Skyrme parameters are also 
determined by the same procedure. The Utot with the 
choice of Vci[A„(^^^Dy)] is almost the same quality com- 
pared to SLy4 (see Table In Fig. [ITl the A.^/A^-^' 
of Sn and Pb isotopes and the Ap/Ap^-* of iV = 50 
and 82 isotones obtained with Skyrme BSklT, LNS and 
SkM* are shown. The 771 and Vo[/^n{^^^^y)] in Table 
mil are used for the Skyrme forces. For comparison, 
the value obtained with the SLy4 force is subtracted; 
5{^r/^r^^){X) = ^r|^''r^\x) - A./A^"^) (SLy4) for 
X = BSkl7, LNS and SkM*. Their (Ami, m^/m) 
are (0.190,0.800), (-0.028,0.780), (-0.164,0.727), and 
(-0.262,0.653) for SLy4, BSkl7, LNS, and SkM* re- 
spectively. The cTtot of the BSkl7 is smallest, and the 
K of SkM* is largest in this work. The LNS parametriza- 
tion was built to match the / dependence of the effec- 
tive masses and the neutron matter EOS predicted by 
Briickner-Hartree-Fock calculation. 

In spite of the variety of Ami and m*, the (5(At-/At^-') 
is small along the isotopic and isotonic chains, except for 
around the subshell closure; iV = 90 for A„, and Z = 38 
and 40 for Ap. This is because the pairing correlations 
are sensitive to the single-particle structure around the 
subshell closure. 

The smah 5(A„/Ai'^^) can be expected due to the 
weak sensitivity of cr„ to rji if the strength Vq is con- 
strained by A^'''^^-' of a specific nucleus [l^l • This is seen 
in Fig. [HI Authors of Ref. [Ill also pointed out that the 
A„ of Sn and Pb isotope chains are insensitive to Ami 
by performing the HFB calculation with the mixed type 
pairing force and various Skyrme forces. On the other 
hand, the fine tuning of 771 is indispensable for the small 
(5(Ap/Ap^^) due to the sensitivity of to r/i. This is 
shown in Fig. [51 and discussed in Ref. [2J] . 

The 772 — 2.5 for SLy4 is commonly used for other 
Skyrme parameters (see Sec. IVII) . This is an approxima- 
tion in our analysis. However, the small 5{/S.t- / as a 
function a means that the difference in the a dependence 
of the pairing gaps due to the different Skyrme forces can 
be small with the fixed 772 — 2.5. 

We refer the pair-DF with the parameters in Table [Till 
as the optimal one for each Skyrme parameterization. 
The optimal pair-DF with pi dependence is constructed 
aiming at unique description of pairing properties toward 
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the neutron drip line. Our prescription is based on the 
phenomenological considerations. In this sense, our con- 
clusion is tentative. However, the optimal pair-DF can 
preserves the good descriptive power of the neutron ex- 
cess dependence of pairing correlations, and provide the 
certain foundation for the further improvement with ex- 
perimental data of nuclei with larger neutron excess. 



ing gaps in finite nuclei based on the phenomenological 
consideration. To obtain the deeper insight to the pi 
terms in the pair-DF, it is interesting to investigate pair- 
ing correlations in asymmetric nuclear matter by com- 
paring the up-to-date calculations with 3-body force and 
correlations beyond the mean-field approximation. This 
analysis is a future subject. 



VIII. CONCLUSION 

We proposed a new pair-DF by introducing the pi de- 
pendence. We emphasized the necessity of both the lin- 
ear and quadratic pi terms in the pair-DF for the global 
description of pairing correlations; namely the depen- 
dence on both the mass number A and the neutron excess 
a={N- Z)/A. 

To optimize the parameters in the pair-DF, we per- 
formed the HFB calculation for 156 nuclei oi A = 118 — 
196 and a < 0.25. By the extensive investigation with 13 
Skyrme parameterizations, we clarify that the pair-DF 
should include the pi dependence in order to take into 
account the effect of the m* and m* in the p-h channel: 
The 771 and m/m* is linearly dependent, and the pairing 
strength Vq linearly increases as a function of m* / m with 
the optimal set of (jyo, ??2)- The relationship between 
the optimal r/i and the splitting of the neutron and pro- 
ton effective masses is also discussed. The Vq is fixed so 
as reproduce the An^'^^ of ^'''^'Dy. With this choice, we 
can obtain the almost minimum value of the total r.m.s. 
deviation between the experimental and calculated pair- 
ing gaps. The different Skyrme forces with the optimal 
pair-DF can give the small difference in the pairing gaps 
toward the neutron drip line. 

In this paper, we concentrated on the analysis of pair- 
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APPENDIX A: CORRELATION ANALYSIS 

We introduce the correlation coefficient r for a data set 
{x, y) = {{xi,yi)} (i = 1, 2, • • • , n). It is defined by 



The coefficient r can take a real value of — 1 < r < 1 . In 
the limit of r = 1 or —1, the data set is linearly depen- 
dent. On the other hand, the correlation between x and 
y is weak if r is close to zero. 
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